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Abstract. Supernova flux and polarization spectra bring vital information on the 
geometry, physical conditions, and composition structure of the ejected matter. For 
some supernovae the circumstellar matter is also probed by the observed spectra. Some 
of this information can be inferred directly from the observed line profiles and fluxes, 
but because of the Doppler broadening and severe line blending, interpretation often 
involves the use of synthetic spectra. The emphasis in this Chapter is on recent results 
obtained with the help of synthetic spectra. 



1 Introduction 

As discussed and illustrated by M. Turatto in Chap. 3, optical spectra serve as 
the basis for the classification of supernovae (SNe). This chapter is concerned 
with extracting information from optical spectra on the geometry, physical con- 
ditions, and composition structure (element composition versus ejection velocity) 
of the ejected matter, and also the circumstellar matter for some SNe. Several 
other chapters of this volume also discuss certain aspects of SN spectroscopy: ul- 
traviolet spectra (Chap. 7 by N. Panagia); the circumstellar and nebular spectra 
of SN 1987A (Chap. 15 by R. McCray); and the hypernova SN 1998bw (Chap. 25 
by K. Iwamoto et al.). 

Some of the information carried by SN spectra can be inferred directly from 
the observed line profiles and fluxes, but because SN ejection velocities are a 
few percent of the speed of light, spectral features generally are blended and 
interpretation often involves comparing observed spectra with synthetic spectra 
calculated for model SNe. The spectrum calculations may be simplified and 
rapid, or physically self-consistent and computationally intensive. Similarly, both 
simple parameterized physical models and numerical hydrodynamical models are 
used. In this chapter we place some emphasis on recent results obtained with 
the help of synthetic spectra. 

In §2 the elements of spectral line formation in SNe are discussed and some 
of the synthetic spectrum codes that are in current use are described. In §§3-6 
recent comparisons of synthetic spectra with observed spectra of SNe of Types 
la, lb, Ic, and II are surveyed. Space limitations prevent us from referring to the 
numerous papers that have been important to the development of this subject; 



many such can be found in 1993 and 1997 review articles 134 ,36]. Here most of 
the references are to, and all of the figures are from, papers that have appeared 
since 1998. A brief discussion of the prospects for the future appears in §7. 
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2 Line Formation and Synthetic Spectrum Codes 

At the time of the supernova explosion, physical scales are relatively small. 
Shortly thereafter the ejected matter expands freely as from a point explosion 
with each matter element having a constant velocity. All structures in the ejecta 
then just scale linearly with t, the time since explosion. The matter is in homol- 
ogous expansion, which has several convenient features: (1) the radial velocity v 
of a matter element is a useful comoving coordinate with actual radial position 
of the element given by r = vt; (2) the density at any comoving point just scales 
as t~ 3 ; (3) the photon redshift between matter elements separated by velocity 
Av, AX = X(Av/c), is time- independent; and (4) the resonance surfaces for line 
emission at a single Doppler-shifted line frequency are just planes perpendicular 
to the observer's line of sight. (Relativistic effects introduce a slight curvature 



If continuous opacity in the line forming region is disregarded, the profile of 
an unblended line can be calculated when the line optical depth T£ (v) and source 
function Sg(v) are specified. Because SN ejection velocities (~ 10,000 km s -1 ) 
are much larger than the random thermal velocities (~ 10 km s _1 ), a photon 
remains in resonance with an atomic transition only within a small resonance 
region. The Sobolev approximation, that the physical conditions other than the 
velocity are uniform within the resonance region, usually is a good one, and it 
allows the optical depth of a line to be simply expressed in terms of the local 
number densities of atoms or ions in the lower and upper levels of the transition: 



where / is the oscillator strength, A M is the line wavelength in microns, td is the 
time since explosion in days, and nt and n u are the populations of the lower and 
upper levels of the transition in cm~ 3 . The term in brackets is the correction for 
stimulated emission. The source function is 
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All of the radial dependence of ti and Si is in the level populations. The specific 
intensity that emerges from a resonance region is 



Spectroscopic evolution can be divided into a photospheric phase when the 
SN is optically thick in the continuum below a photospheric velocity, and a sub- 
sequent nebular phase during which the whole SN is optically thin in the contin- 
uum. In the photospheric phase line formation occurs above the photosphere and 
in the nebular phase throughout the ejecta. There is, of course, no sharp division 
between the two phases. However, spectral synthesis modeling techniques in the 
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photospheric and nebular limits can make use of different approximations which 
are adequate for those limits. 

In the remainder of this section we first discuss line formation and synthetic 
spectrum codes in the photospheric and nebular phases under the assumptions 
of spherical symmetry and negligible circumstellar interaction. We then consider 
the effects of circumstellar interaction, and SN asymmetry and polarization spec- 
tra. 



2.1 The Photospheric Phase 

The elements of spectrum formation in the photospheric phase have been dis- 
cussed and illustrated at length elsewhere |66|1 . Here we will only briefly summa- 
rize those elements and then discuss spectrum synthesis codes. 

During the photospheric phase a continuum radiation field is emitted by a 
photosphere which can be idealized as an infinitely thin layer. Above the photo- 
sphere the radiation interaction with continuous opacity is small. Line opacity 
on the other hand can be very large for the strongest lines. The large Dopplcr 
shifts spread line opacity over a large wavelength interval increasing the effect 
of strong lines compared to a static atmosphere where such strong lines saturate 
and can only affect radiation in a narrow wavelength interval. The cumulative 
effect of many lines, strong and weak, can create a quasi-continuous opacity in 
the Eulerian frame. This effect has been called the "expansion opacity" |6£,107|; 



it dominates in the ultraviolet, where it effectively pushes the photosphere out 
to a larger radius than in the optical. Full, co-moving, line-blanketed calcula- 
tions automatically account for this effect ||; introducing the observer's frame 
expansion opacity into comoving frame calculations is incorrect. 

In the optical, which is usually the chief focus of analysis, the spectrum is 
characterized by P Cygni lines superimposed on the photospheric continuum. 
The P Cygni profile has an emission peak near the rest wavelength of the line 
and a blueshifted absorption feature. The peak may be formed in part by true 
emission or by line scattering into the line of sight of photons emitted by the 
photosphere. The emission peak would tend to be symmetrical about the line 
center wavelength if not for the blueshifted absorption. The absorption is formed 
by scattering out of the line of sight of photospheric photons emitted toward 
the observer. Since this occurs in front of the photosphere, the absorption is 
blueshifted. At early times the ejecta density is high, the photosphere is at high 
velocity, and the line opacity is strong out to still higher velocities. As expansion 
proceeds, the photosphere and the region of line formation recede deeper into the 
ejecta. The P Cygni line profile width thus decreases with time. The minimum 
of the absorption feature of weak lines tends to form near the photospheric 
velocity, thus weak lines (e.g., weak Fe II lines) can be used to determine the 
photospheric velocity's time evolution. The recession of the photosphere exposes 
the inner ejecta and permits its analysis. 

In understanding P Cygni line formation in SNe it is conceptually useful, 
and often in simple parameterized radiative transfer modeling a good physical 
approximation, to consider the line-photon interaction as being pure resonance 
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scattering; i.e., no photon creation or destruction above the photosphere. Then 
the source function of an unblended line is just Se(v) — W(v)I p )x t, where W(v) 
is the usual geometrical dilution factor (9{| and / p hot is the continuum specific 
intensity (assumed angle-independent) radiated by the photosphere. Given the 
pure resonance scattering approximation for an unblended P Cygni line, the 
emission is formed just by scattering into the line of sight and the absorption 
just by scattering out of the line of sight. 

Unfortunately, in general there is strong line blending. A photon that is scat- 
tered by one transition can redshift into resonance with another, so the influence 
of each transition on others of longer wavelength must be taken into account. 
This multiple scattering corresponds to the observer's line blending. Calcula- 
tions show that in blends, absorptions trump emissions: i.e., when resolved, the 
absorption minima of blends tend to be blueshifted by their usual amounts, but 
the emission peaks do not necessarily correspond to the rest wavelengths of the 
lines of the blend. Thus absorption minima usually are more useful than emission 
peaks for making line identifications during the photospheric phase. 

If one has made the resonance scattering approximation for simple modeling, 
it turns out that the effect of line blending on source functions for some lines is 
often effectively very small. From a conceptual point of view it is convenient to 
think of Se(v) = VF(u)^ P hot even for strongly blended lines. In general the effect 
of line blending on the emergent spectrum is significant and cannot be neglected. 

A special case of a P Cygni line that has become of interest is a "detached" 
line: i.e, a line that has a significant optical depth only above some detachment 
velocity that exceeds the velocity at the photosphere. A detached line consists 
of a flat inconspicuous emission peak and an absorption having a sharp red edge 
at the blueshift corresponding to the detachment velocity (Fig 1). We shall later 
refer to the concept of detached lines in discussing particular observed spectra. 

It is possible to extract some information from spectra by simple direct analy- 
sis. Recently a method for inverting an unblended P Cygni line profile to extract 
the radial dependences of the line optical depth and the source function has 
been developed |^|. But the complex blending of P Cygni lines (which often 
occurs), and the need to account for effects that cannot treated by direct ap- 
proaches make analysis by synthetic spectrum modeling essential. Rather than 
discuss modeling in general we will briefly describe three synthetic spectrum 
codes that are frequently used, proceeding from simple to complex. Example 
synthetic spectra from these codes will appear in §§3 - 6. 

The simplest code is the fast, parameterized SYNOW code. Technical details of 
the current version are in [p8| . The basic assumptions are: spherical symmetry; a 
sharp photosphere that emits a blackbody continuous spectrum; and the Sobolev 
approximation with a resonance scattering source function. SYNOW does not do 
continuum transport, it does not solve rate equations, and it does not calculate or 
assume ionization ratios. Its main function is to take line multiple scattering into 
account so that it can be used in an empirical spirit to make line identifications 
and determine the velocity at the photosphere and the velocity intervals within 
which the presence of each ion is detected. For each ion that is introduced, 
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Fig. 1. A SYNOW synthetic spectrum (dotted line) that has v p hot = 10,000 km s" 1 and 
hydrogen lines detached at 20,000 km s" 1 where r(Ha)=2, is compared with a synthetic 
spectrum (solid line) that has v v hot = 10, 000 km s~ and undetached hydrogen lines 
with r(Ha) = 10 at the photosphere. (From Jig}]) 



the optical depth at the photosphere of a reference line is a free parameter, 
and the optical depths of the other lines of the ion are calculated assuming 
Boltzmann excitation. Reference lines are generally chosen as the strongest line 
in the optical for a particular species. Line optical depths are assumed to decrease 
radially following a power law or an exponential. When deciding which ions to 
introduce, use is made of the results of [5(| who presented plots of LTE Sobolev 
line optical depths versus temperature for six different compositions that might 
be expected to be encountered in SNe, and sample SYNOW optical spectra for 45 
individual ions that are candidates for producing identifiable spectral features 
in SNe. When fitting to an observed spectrum, the important parameters are 
the reference line optical depths, the velocity at the photosphere, and whatever 
maximum and minimum (detachment) velocities may be imposed on each ion. 

Somewhat more complex synthetic spectrum calculations are carried out with 
what we will refer to as the ML MONTE CARLO code (ML for Mazzali and Lucy), 
technical details of which are in |)l| and references therein. This code also as- 
sumes a sharp photosphere, makes the Sobolev approximation, and does not 
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solve rate equations. The main differences from SYNDW are: (1) an approximate 
radiative equilibrium temperature distribution and internally consistent ioniza- 
tion ratios are calculated for an assumed composition structure; (2) electron 
scattering in the atmosphere is taken into account; and (3) the original assump- 
tion of resonance scattering has been replaced by an approximate treatment of 
photon branching, to allow an absorption in one transition to be followed by 
emission in another pBfl . Several circumstances in which the effects of photon 
branching are significant are discussed in (9^] . 

Very detailed calculations are made with the multi-purpose synthetic spec- 
trum and model atmosphere code called PHOENIX. Technical details of PHOENIX 
are in fi3"[ (see also f| ) . The basic assumptions are spherical symmetry and time 
independence. The Sobolev approximation is dispensed with, and calculations 
are carried out in the comoving frame with all special relativistic effects. The 
aim of PHOENIX is to take all of the relevant physics into account as fully as 
possible. PHOENIX treats continuum transport, solves the NLTE rate equations 
for a large number of ions and a very large number of atomic levels with gamma 
ray deposition and nonthermal excitation taken into account, and determines a 
self-consistent radiative equilibrium temperature structure. The lower boundary 
condition can be chosen to be either diffusive or nebular. 

One may ask why simpler codes are needed when PHOENIX-like codes can 
be implemented. PHOENIX-like codes demand many hours of computation and 
are not optimal for explorations of parameter space or for gaining insight into 
how simple features of modeling affect outcomes. The simpler codes are tools for 
rapid exploration and gaining insight. Such codes, however, cannot give defini- 
tive determinations because of their simplified physics. PHOENIX-like calculations 
must be the basis for ultimate decisions about the viability of SN models. 



2.2 The Nebular Phase 

Our discussion of line formation in the nebular phase will be brief since the 
subject is discussed in Chap. 15, in the context of SN 1987A. (See also [pi].) 

In the absence of circumstellar interaction, the nebular phase is powered by 
radioactive decay, primarily 56 Co (77.23 day half- life), the daughter of explosion- 
synthesized 56 Ni (6.075 day half-life). At very late times longer lived radioactive 
species such as 57 Co and 44 Ti become important. 

The decay gamma rays deposit their energy by Compton scattering off elec- 
trons to produce nonthermal fast electrons whose energy quickly goes into atomic 



ionization and excitation, and Coulomb heating of the thermal electrons [ 12C ] . 
Particularly for low mass SNe (Types Ia,b,c) the unscattered escape of gamma rays 
(which increases as the ejecta thins) causes the kinetic energy of positrons from 
the 56 Co decays to become an important energy source: about 3% of the 56 Co 
decay energy is in the form of positron kinetic energy. The positrons are much 



more strongly trapped than the gamma rays — but not completely trapped [101]. 
The positron kinetic energy goes into creating fast electrons leading to the same 
deposition processes as the gamma rays. Optical emission lines are formed by 
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recombination, collisional excitation, and fluorescence. P Cygni profiles of per- 
mitted lines that form by scattering in the outer layers may be superimposed on 
the emission line spectrum. 

Because of low continuous opacity and increasingly low line opacities, pho- 
tons emitted by true emission processes are increasingly unlikely to be scattered 
again. Thus non-local radiative transfer effects in the nebular phase tend to be- 
come negligible. (Local trapping in an optical thick line can be treated easily 
using a Sobolev escape probability.) Thus (absent circumstellar interaction) it 
is somewhat more feasible than in the photospheric phase to extract informa- 
tion without making synthetic spectrum calculations. This was especially true 
of SN 1987A, which because of its proximity was well observed in the optical 
and the infrared long into the nebular phase when the characteristic velocity of 
the line forming region was only ~ 3000 km s^ 1 and blending was not too severe 
(see Chap. 15). Most other non-circumstellar interacting SNe are observed only 
for a year or so after explosion and line blending is more of a problem. 

Without a photosphere, there is no continuum flux to be scattered out of the 
line of sight even if there are optically thick lines. Hence the characteristic line 
profile of the nebular phase is not a P Cygni profile, but an emission line that 
peaks near the rest wavelength and for several reasons may have an extended red 
wing |^3|,|o) . If the line has significant optical depth down to low velocity, the 
emission line has a rounded peak. If instead the line forms in a shell, the emission 
has a boxy (flat top) shape with the half width of the box corresponding to the 
minimum velocity of the shell. 

As for the photospheric phase, nebular-phase synthetic spectrum calculations 
of various levels of complexity prove to be useful. Relatively simple parameter- 
ized approaches that can be used as diagnostic tools are described by JlC|fL| . 
A one-zone code that was developed in its original form by |11C[ and will be 
referred to here as the RL NEBULAR code (RL for Ruiz-Lapuente and Lucy) is 
frequently used (e.g. |)2,D3 109 ). The input parameters are the time since ex- 



plosion and the mass, composition, and outer velocity of the nebula. Heating is 
calculated from the 56 Co decay rate, and cooling is by radiative emission. The 
electron density, temperature, NLTE level populations and emission line fluxes 
are determined simultaneously by enforcing statistical and thermal equilibrium. 
Codes for calculating synthetic spectra of hydrodynamical models are described 
by |jo| and ||| and references therein. These multi-zone codes calculate the 
nonthermal heating, ionization equilibrium, and the NLTE level populations for 
a large number of ato mic l evels. Synthetic spectra of the CO molecule have also 
been calculated |3^ , ^5| , pT8| . For a description of additional nebular phase codes 
that have been used for SN 1987A, see Chap. 15. 

The photospheric phase forbids a direct view into the inner regions of the 
ejecta. The nebular phase permits this direct view, but, on the other hand, the 
outer layers tend toward invisibility. The photospheric phase demands a powerful 
radiative transfer technique. The advanced nebular phase requires almost none 
at all, at least in the optical; the ionizing photons in the UV may need a detailed 
radiative transfer treatment and, of course, the gamma ray and positron energy 
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depositions need a non-local transfer treatment probably best done by Monte 
Carlo. Both phases require NLTE treatments and extensive atomic data. The 
two phases offer complementary insights and challenges for analysis. 



2.3 Circumstellar Interaction 

The physics of the radiative and hydrodynamical interactions between SNe 
and their circumstellar matter (circumstellar interaction: CS1) are discussed in 
Chap. 14 by R.A. Chevalier & C. Fransson. No convincing detection of SN la 
CSI has yet been made in any wavelength band. CSI also does not seem to affect 
the optical spectra of typical SNe lb, SNe Ic (but for an exception see and 
many SNe II. Some SNe II and perhaps some hyperenergetic SNe Ic are affected, 
and in extreme cases dominated, by CSI effects. 

In the idealized case, the wind of the SN progenitor star has a constant ve- 
locity and a constant mass loss rate, and therefore an r~ 2 density distribution. 
After the explosion, the circumstellar matter can be heated and accelerated by 
the photons emitted by the SN. Subsequent hydrodynamic interaction between 
high velocity SN ejecta and low velocity circumstellar matter generates photon 
emission that further affects the velocity distribution and the ionization and exci- 
tation of both the circumstellar matter and the SN ejecta. These interactions can 
produce a rich array of effects on the optical spectra. During the photospheric 
phase narrow emission, absorption, and P Cygni lines formed in lo w velocit y cir- 
cumstellar matter may be superimposed on the SN spectrum (|^, 112, 116 1. The 



SN resonance scattering features may be "muted" by the external illumination 
of the SN line forming region [M . In extreme cases the circumstellar matter may 
be optically thick in the continuum so that the spectrum consists only of circum- 
stellar features, perhaps including broad wings produced by multiple scattering 
off circumstellar electrons ]2l]| . During the nebular phase CSI can power boxy 
emission from high velocity SN matter ]42| and at very late times, even decades, 
CSI can allow the detection of SNe that would otherwise be unobservably faint 

ID- 

Because the simple r = vt velocity law does not apply to circumstellar mat- 
ter, and blending of circumstellar lines generally is not as severe as it is for 
SN lines (and spherical symmetry often is not an adequate approximation for 
circumstellar matter), few calculations of synthetic spectra for CSI have been 
carried out so far. Line profiles for various velocity and density distributions have 
been calculated flippy and some PHOENIX calculations for a constant velocity 
circumstellar shell have been carried out J7(|. Most of the analysis of optical CSI 
features has involved extracting information on velocities directly from individ- 
ual line profiles, and on physical conditions and abundances from emission line 
fluxes |T7|J2l,|32lM. 



2.4 Asymmetry and Polarization Spectra 



It has been the heuristic and hopeful assumption that SNe are essentially spher- 
ically symmetric. But from SN polarimetry, especially spectropolarimetry, it is 
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now clear that all types of SNe other than SNe la usually do exhibit some kind 
of significant asymmetry (e.g, [61,78 130, 132| ). 

The polarization of SN flux in most cases arises from their atmospheres where 
continuous opacity is dominated by electron scattering in the optical and near 
infrared. Electron scattering is polarizing with polarization (for incident unpo- 
larized radiation) varying from 100 % for 90° scattering to % for forward and 
backward scattering according to the Rayleigh phase matrix (l6| . The polariza- 
tion position angle is perpendicular to the scattering plane defined by incident 
and scattered light. The emergent flux from a deep electron scattering atmo- 
sphere (either plane-parallel or geometrically extended) tends to exhibit a posi- 
tion angle aligned perpendicular to the normal to the surface and polarization 
increasing monotonically from % at an emission angle of 0° (to the normal) to a 
maximum at an emission angle of 90° . For the plane-parallel case the maximum 
is only 11.7% fL6f , but for a highly extended atmosphere it can approach 100 % 
[ fi"5| from the limb, if it can be resolved. SN atmospheres fall into the extended 
class, and so beams from their limb should be highly polarized though mostly not 
at 100 %. Beams from the photodisk region (the projection of the photosphere on 
the sky) should be much less polarized. If SNe were spherically symmetric, the 
integrated flux (which is all that can be observed from extragalactic SNe) would 
have zero polarization since all polarization alignments would contribute equally 
and thus cancel producing zero net polarization. Because of cancellation, the 
predicted net continuum polarizations are low (a few percent: see below) even 
for quite strong asymmetries. 

Line scattering is usually depolarizing p2[ |. The simplest line polarization 
profile is an inverted P Cygni profile that is common although not universal 
in the SN 1987A data (63) and is seen in some of the later spectropolarimetry 
[ fL30| . In this profile a polarization maximum is associated with the P Cygni 
absorption feature of the flux spectrum and a polarization minimum is associated 
with the P Cygni flux emission feature. The inverted P Cygni profile has a 
natural explanation. The polarized flux at any wavelength tends to come from 
90° electron scattering. At the flux absorption, unpolarized flux directly from 
the photosphere is scattered out of the line of sight by a line: thus, the electron 
polarized flux is less diluted and a polarization maximum arises. At the flux 
emission, the line scatters into the line of sight extra unpolarized flux diluting 
the polarization due to the electron scattered flux. Almost any shape asymmetry 
should tend to give the inverted polarization P Cygni profile. Thus, the inverted 
P Cygni profile is somewhat limited as a direct analysis tool for asymmetry. 

Although one must be cautious in generalizing because of the paucity and 
inhomogeneity of published data, it seems that most core-collapse SNe that have 
been observed well enough to detect polarization of order a few tenths of per- 
cent (of order 40 to date) have shown at least that much continuum polarization. 
Some of this can be interstellar polarization (ISP) due to intervening dust, but 
intrinsic polarization seems to be present in many cases. The intrinsic polariza- 
tion reveals itself through time dependence and, in spectropolarimetry, by line 
polarization features or by deviation from the ISP wavelength dependence [114]. 
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Unfortunately, correcting for parent-galaxy ISP can only be done from an anal- 
ysis of the SN data itself because there is no other bright nearby point source. 
Even the Galactic ISP correction may often have to be determined using the 
SN data since ISP can vary rapidly with position and depth and there may be 
no suitable Galactic star near the angular position of the SN. Techniques for 
correcting for ISP have been developed |6l]|33|,[77j . 

Fig. 2 shows a sample of spectropolarimetry from a range of SN types overlap- 
ping the corresponding relative flux spectra. The polarization has been corrected 
for ISP, but in most cases only with the intention of highlighting the line po- 
larization features: the overall polarization level is not trustworthy. Some of the 
line polarization features are consistent with the inverted P Cygni profile. How- 
ever, line polarization features are badly distorted by ISP, and thus reliable ISP 
corrections are needed to determine the true line polarization profiles as well as 
the true continuum polarization. 

For core-collapse SNe the intrinsic continuum polarization varies from prob- 
ably ~ 0.0 % up to perhaps ~ 4 % |7^ , |l30| : a representative value would be 1 %. 
There may be a class of core-collapse SNe (those with massive hydrogen en- 
velopes) that have nearly zero polarization at early times [Q , and so are highly 
spherical. Apparent trends are that polarization increases with decreasing en- 
velope mass (i.e., SNe II-P least polarized, SNe Ic most), and that polarization 
increases with time until the electron optical depth becomes small |33]]76,13C]. 
The position angle of polarization tends to be nearly constant in time and wave- 
length [6^,76 130 , which is naturally accounted for by axisymmetry in the ejecta. 

For SNe la most observations (of order 15 to date) have only been able to 
place an upper limit of about 0.3 % continuum polarization (not distinguishing 
intrinsic from ISP in most cases) |132|. One normal event, SN 1996X, has a 



marginal detection of intrinsic polarization of ~ 0.3 % [129]. Spectropolarimetry 
of SN 1997dt also may show the signatures of intrinsic polarization (Fig. 2). So far 
it seems that normal SNe la are not very polarized and are quite spherically sym- 
metric: this is consistent with the observational homogeneity that makes them 
such useful tools for cosmology. On the other hand, the subluminous SN 1991bg- 
like event SN 1999by had intrinsic continuum polarization reaching up to perhaps 
0.8% [pl|| . Based on a sample of one, subluminous SNe la may be significantly 
asymmetric. 

Various models of asymmetry have been presented to interpret the polariza- 
tion data. The conventiona l mo del is axisymmetric ellipsoidal asymmetry, either 
prolate or oblate |^6| , |62|p^ , |ll5|| . Even when viewed equator-on, considerable el- 
lipsoidal asymmetry is required to account for SN-like continuum polarization. 
For example, Monte Carlo calculations suggest that iso-density contour axis ra- 
tios of or der 1 .2:1 and 2.5:1 are needed for polarizations of 1 % and 3 %, respec- 
tively [58 130 . Late -time imaging of SN 1987A shows asymmetries suggesting 
axis ratios of 2:1 |p7 131 1. Thus large organized asymmetries something like the 
ellipsoidal models may account for most SN polarization. Other suggested mod- 
els of polarizing large-scale asymmetry are asymmetric clumpy ionization ]lg| ], 
or scattering of SN light off a nearby dust cloud 1 12S ] or off bipolar jets |55| . 
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Fig. 2. Spectropolarimetry (thin lines) of various SN types, with the corresponding 
flux spectra (thick lines), all obtained within two months of explosion. Figure courtesy 
of D.C. Leonard 

Spectropolarimetry plays an extremely important role in polarimetric obser- 
vations for two reasons: (1) it is very useful for establishing intrinsic polarization 
and correcting for ISP, and (2) it is probably vital for understanding the SN 
asymmetry. Currently the most advanced synthetic polarization spectrum anal- 
yses are done with a 3-dimensional (3D) NLTE code that uses Monte Carlo 
radiative transfer in the outer region where emergent polarization is formed 
| j57[|6l| . The NLTE component of these calculations is necessarily more limited 
than in ID calculations since 3D radiative transfer calculations are computa- 
tionally very demanding. The ellipsoidal asymmetry used in the calculations of 
[ |6l| is parameterized and then fitted to the flux and polarization observations. 
Explosion model asymmetries can also be used |57| . Because of the large param- 
eter space to explore in investigating SN asymmetry, less elaborate codes than 
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that of |37j and will also continue to be useful. Such codes will mostly rely 
on Monte Carlo radiative transfer because of its flexibility and robustness. 

In addition to the shape asymmetries discussed above, the ejected matter 
also may be "clumped" (see Chap. 15). Both observation and theory indicate 
that macroscopic (not microscopic) mixing is ubiquitous in SNe. Clumping is 
difficult to detect via polarization (although see p8|), but it can affect both the 
photospheric [||§3g| and nebular [|l0y24y3^,|7^,g8|]^| flux spectra. A version of 
SYNOW for calculating flux spectra with clum ps (CLUMPYSYN) has been developed 
for analysis of the photospheric phase [122]. 

Circumstellar matter may and probably often does, have both global asym- 
metry and clumping jgJQ^Q . 



3 SNe la 

SNe la are thought to be thermonuclear disruptions of accreting or merging 
white dwarfs. They can be separated on the basis of their photospheric spectra 
into those that are normal and th ose t hat are peculiar. Normals such as the 
recently well observed SNe 1996X (m) and 1998bu ||,|7) have P Cygni fea- 



tures due to Si II, Ca II, S II, O I, and Mg II around the time of maximum 
brightness, and develop strong Fc II features soon after maximum. Peculiar, 



weak SN 1991bg-like events such as SNe 1998de |03|, 1997cn JL25J, and 1999by 



|f4[|61 ,127| have, in addition, conspicuous low excitation Ti II features. Peculiar 
powerful SN 1991T-like events such as SNe 1997br Q and 2000cx @ have 
conspicuous high excitation Fe III features at maximum and the usual SN la 
features develop later. Events such as SN 1999aa appear to link the 1991T- 
like events to the normals pj] |. Nebular spectra of SNe la are dominated by 
collisionally excited forbidden emission lines of [Fe III], [Fe II], and [Co III]. 
SN la ejecta basically consist of a core of iron peak isotopes (initially mostly 
56 Ni) surrounded by a layer of intermediate-mass elements such as silicon, sul- 
fur, and calcium. Inferring the composition of the outermost layers is difficult 



] p9[|49| , |92| , |95[ because the lines are broad and line formation occurs in these layers 
for only a short time after explosion. 

Current viable hydrodynamical models for SNe la include deflagrations and 
delayed detonations (see Chap. 12 by S.E. Woo sley). The composition structure 
of the classic deflagration model W7 [ 104 |l21 ] can account reasonably well for 



most of the photospheric and nebular spectral features of normal SNe la (Figs. 
3 and 4). On the other hand, near infrared spectra appear to contain valuable 
diagnostics of the composition structure |jl^,|3J9|] and may provide evidence in 



favor of delayed detonation models [133] 



SNe la can be arranged in a spectroscopic sequence ranging from the powerful 
SN 1991T-like events through the normals to the weak SN 1991bg-like events. 



The sequence usually is based on the R(Si II) parameter [ 105 1 - the ratio of the 
depth of an absorption feature near 5800 A to that of the deep Si II absorption 
feature near 6100 A. The ratio increases with decreasing temperature owing to 
the increasing contribution of numerous weak Ti II lines to the 5800 A absorption 
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Fig. 3. The spectrum of the Type la SN 1994D three days before maximum brightness 
(solid line) and a PHOENIX synthetic spectrum for the deflagration model W7 17 days 
after explosion (dashed line) (From [[72| ) , 



[ p4| (Fig. 5). The velocity at the outer boundary of the iron peak core, as inferred 
from nebular spectra, also varies along the spectral sequence |)0| . The principle 
physical variable along the sequence is thought to be the mass of ejected 56 Ni in 
a constant total ejected mass which is near the Chandrasekhar mass. 

The lack of a clear correlation between R(Si II) and the blueshift of the 6100 A 
feature shows that the diversity among SN la spectra actually is at least two 
dimensional 51 . Blueshift differences among events that have similar values of 
R(Si II) appear to be caused by differences in the amount of mass that is ejected 
at high velocity, - 15, 000-20, 000 km s" 1 JnJ (Fig 6). One possibility is that the 
high velocity events are delayed detonations while the others are deflagrations. 

The influence of metallicity on SN la spectra is of special interest in connec- 
tion with the use of high redshift SNe la as distance indicators for cosmology. 
The main effects are in the ultraviolet but there also are some mild effects in the 
optical @Q. 

Deflagrations are expected to produce macroscopic mixing and clumping in 
the ejecta. The first application of the CLUMPYSYN code has shown that the 
uniformity of the observed depth of the 6100 A Si II absorption in SNe la near 
maximum brightness provides limits on the sizes and numbers of clumps [122]. 

Although the polarization of most SNe la is low, the significant polarization 
of the SN 1991bg-like event SN 1999by (Fig. 7) raises the question of whether the 
production of weak SNe la may involve rapid white dwarf rotation or merging 

iBil. 




Fig. 4. The spectrum of the Type la SN 1996X 298 days after maximum brightness 
(solid line) is compared with a synthetic spectrum calculated with the RL NEBULAR 
code (dotted line) (From 1 1 1 1 ] ) 



4 SNe lb 

SNe lb are thought to result from core collapse in massive stars that have lost all 

or almost all of their hydrogen envelopes. The photospheric spectra contain the 

usual low excitation SN features such as Ca II and Fe II together with strong 

He I lines that are nonthermally excited by the decay products of 56 Ni and 
56 r 



Co pq,|84,119|. Recent work on photospheric spectrum calculations includes a 



comparison of one NLTE synthetic spectrum of a hydrodynamic model with a 



spectrum of SN 1984L |13fJ, a SYNOW study of line identifications in SN 1999dn 



[ p0| , and a comparative SYNOW study of the spectra of a dozen SNe lb 13 . 

SYNOW spectra can match the observed spectra rather well (Fig. 8). Consistent 
with the observational homogeneity of SNe lb found by (8^] , the sample studied 
by |l3]| obeys a surprisingly tight relation between the velocity at the photosphere 
as inferred from Fe II lines and the time relative to maximum brightness. The 
masses and kinetic energies of the events in the sample appear to be similar, and 
not much room is left for any influence of departures from spherical symmetry on 
the velocity at the photosphere. After maximum brightness the He I lines usually 
are detached, but the minimum velocity of the ejected helium is at least as low 
as 7000 km s" 1 (Fig. 9). The spectra of SNe 2000H, 1999di, and 1954A contain 
detached hydrogen absorption features forming at ~ 12, 000 km s _1 (Fig. 9), and 
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Fig. 5. SYNOW synthetic spectra with only Si II lines {heavy line) and also with Ti II 
lines of various strengths (thin lines). The 580 nm absorption is strongly affected by 
Ti II lines (From G3). 



hydrogen appears to be present in SNe lb in general, although in most events 
it becomes too weak to identify soon after maximum brightness. The hydrogen 
line optical depths used to fit the spectra of SNe 2000H, 1999di, and 1954A are 
not high, so only a mild reduction would be required to make these events look 
like typical SNe lb. Similarly, the He I line optical depths in typical SNe lb are 
not very high, so a moderate reduction would make them look like SNe Ic. 

The number of SNe lb for which good spectral coverage is available is still 
small. More events need to be observed to explore the extent of the spectral 
homogeneity and to determine whether there is a continuum of hydrogen line 
strengths. Also needed are detailed NLTE calculations for hydrodynamical SN lb 
models having radially stratified composition structures, to determine the hydro- 
gen and helium masses and the distribution of the 56 Ni that excites the helium. 

To our knowledge there have not been any synthetic spectrum calculations 
for the SN lb nebular phase since 1989 Jig] . Most of the recent effort has been 
devoted to inferring the composition from emission line profiles and fluxes p9| . 
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Fig. 6. PHOENIX synthetic spectra for the delayed detonation model CS15DD3 (dashed 
lines) and observed spectra of the high- velocity Type la SN 1984A (solid lines). (From 
0) 



The oxygen mass, evidence for hydrogen, and evidenc e for asymmetry from the 
nebular spectra of SN 1996N have been discussed by [117]. 



5 SNe Ic 



SNe Ic are thought to result from core collapse in massive stars that either 
have lost their helium layer or fail to nonthermally excite their helium. If the 
helium layer is gone, the outer layers of the ejected matter are expected to 
be mainly carbon and oxygen. SNe Ic are spectroscopically more diverse than 
SNe lb |28|,^9j . The photospheric spectra are dominated by the low excitation SN 
features that appear in SNe lb. The spectra of the best studied or dina ry Type Ic, 
SN 19941, can be reasonably well matched by SYN0W spectra flOCl (Fig. 10). 
Discrepancies between the SN 19941 spectra and PHOENIX spectra indicate that 
the appropriate hydrodynamical model for SN 19941 has not yet been found ||. 
The important issue of whether helium is detectable in SN Ic spectra is not yet 
resolved, because an absorption near 10,000 A that sometimes is attributed to 
H e I A 10,830 could be produced by one of several other ions (see Fig. 10 and 
[5 10C[), and identifications of very weak optical lines j?6) are difficult to confirm 
|9|. 

The peculiar, hyperenergetic Type Ic SN 1998bw that was associated with 
a gamma ray burst is discussed in Chap. 25. SN 1997ef was a less extreme 
example of a hyperenergetic SN Ic. The spectra of these events, although diffi- 
cult to interpret because of the severe Doppler broadening and blending, can be 
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Fig. 7. Flux spectra (top panel, black line) and raw and smoothed polarization spectra 
(second panel, light blue and black lines) of the peculiar weak Type la SN 1999by near 
maximum brightness, are compared with synthetic spectra computed for an oblate 
ellipsoidal model with axis ratio 1.17 and three inclinations of the symmetry axis from 
the line of sight. (From |pH , where further explanation of this figure can be found) 



matched fairly well by the SYNDW Q and ML MONTE CARLO codes Q (Fig. 11), 
and the resulting spectroscopic estimates of the kinetic energy (in the sph erical 
approximation) are far in excess of the canonical 10 51 ergs 12 9^ , |96| , ^06| . De- 
tailed NLTE synthetic spectrum calculations have not yet been carried out for 
the photospheric spectra of hyperenergetic SNe Ic, and the only recent synthetic 
spectrum calculations for nebular phase SNe Ic have been for SN 1998bw with 
the RL NEBULAR code ||| (see Chap. 25). 
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Fig. 8. The spectrum of the Type lb SN 1999dn 17 days after maximum brightness 
(solid line) is compared with a SYNOW synthetic spectrum (dotted line) that has v p hot = 
6000 km s" 1 . (From (l3)) 

The Type Ic SN 1999as, one of the most luminous reliably measured SNe 
to date, contained highly blueshifted (~ 11,000 km s _1 ) but narrow (~ 2000 
km s _1 ) absorption features ]52]| that cannot be produced by spherically sym- 
metric SN ejecta. Whether these unusual features are produced by an ejecta 
clump in front of the photosphere or by circumstcllar matter accelerated by the 
SN is not yet clear. 

Detailed NLTE calculations for hydrodynamical models of typical and hy- 
perenergetic SNe Ic are needed, to constrain the helium mass and the ejected 
mass and kinetic energy. 



6 SNe II 

SNe II are the most spectroscopically diverse of the SN types, partly because 
the optical effects of CSI range from negligible in some events to dominant 
in others. In the absence of CSI and the presence of a substantial hydrogen 
envelope (SNe II-P and II-L) the spectrum evolves from an almost featureless 
continuum when the temperature is high to one that contains first hydrogen 
and helium lines and then also the usual low excitation lines as the temperature 
falls. The photospheric spectra of SN 1987A were, of course, extensively studied 



[ pl|j55|j66 ,113|. Recently it has been shown that the observed profile of the He I 
A10830 line in SN 1987A is a sensitive probe of the degree of mixing of 56 Ni 




Fig. 9. For a sample of 12 SNe lb, the minimum velocity of the He I lines (filled squares 
when undetached, filled diamonds when detached) and the minimum velocity of the 
hydrogen lines (open circles, always detached) are plotted against time after maximum 
brightness. The heavy curve is a power law fit to the velocity at the photosphere as 
determined by Fe II lines. (From [1^]) 



out to relatively high velocities, ~ 3000 km s 1 133] (see also p4j for similar 



work on the Type II-P SN 1995V). PHOENIX calculations |102| show that the 
hydrogen lines of SN 1987A also are a sensitive probe of the nickel mixing and 
may require some 56 Ni as fast as 5000 km s _1 . As discussed in these three papers, 
quantitative conclusions are sensitive to the details of the macroscopic mixing 
of the hydrogen, helium, and nickel. 

SN 1999em is a recent SN II-P that has been extremely well observed and 
that is not complicated by CSI or large departures from spherical symmetry 
[[76). HST ultraviolet spectra @] of this event are discussed in Chap. 7. Two 
extensive independent sets of optical spectra and photometry have been used to 
apply the expanding photosphere method for determining the distance |4^j75|] . 
SYN0W and PHOENIX synthetic spectra have been compared to the photospheric 
spectra [||. Complicated observed P Cygni profiles of hydrogen and He I lines 
are well reproduced by the PHOENIX calculations when the density gradient is 
sufficiently shallow (Fig. 12). The total extinction can be inferred by comparing 
the observed and synthetic spectra at early times when the temperature is high 
§■ 

The unusually weak Type II-P SN 1997D has attracted considerable attention 
in connection with the possibility of black hole formation and fallback. The ML 
MONTE CARLO code has been used to fit the spectrum with an extremely low 



velocity at the photosphere of only 970 km s 1 [123]. It has been pointed out that 



because of the unusually low ionization, Rayleigh scattering may have significant 
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Fig. 10. The spectrum of the typical Type Ic SN 19941 7 days after maximum bright- 
ness is compared with three SYNOW synthetic spectra that have v p hot = 10, 000 km s _1 . 
To account for the observed infrared feature, the synthetic spectra include He I lines 
detached at 15,000 km s _1 (top); undetached C I lines and He I lines detached at 1 8,000 
km s^ 1 (middle); and Si I lines detached at 14,000 km s (bottom). (From [100]) 



effects on the photospheric spectrum [£5[ . The low expansion velocity offered a 
unique opportunity to identify lines in the nebular spectra [Q . Nebular synthetic 
spectra have been calculated (29] and it has been pointed out that the Sobolev 
approximation may break down in the early nebular phase because Ha damping 
wings can have significant effects |2(| . 

SNe lib contain conspicuous hydrogen lines at early times but not at late 
times, indicating that these events lost nearly all of their hydrogen before ex- 
ploding. In the well studied SN 1993J, He I lines became conspicuous after 
maximum brightness, as the hydrogen lines fade d. The photospheric and nebu- 
lar spectra were modeled in detail (J2,5S,30 126 and references therein). More 
recent work on SN 1993J has concentrated on inferring properties of the clumps 
and the nature of the CSI directly from line profiles and fluxes . A SYNOW 
l ine- identification study of extensive observations of the Type lib SN 1996cb 
[ |108| has been carried out [p9| . 

SNe Iln have narrow lines that form in low velocity circumstellar matter. Op- 
tical and infrared spectra of SN 1998S have been intensively studied ]l|, |32] , [70|j77| , p3] | 
HST ultraviolet spectra jTy] are discussed in Chap. 7. PHOENIX synthetic spectra 
calculated for a simple model of the circumstellar shell accounted reasonably 
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Fig. 11. The spectrum of the hyperenergetic Type Ic SN 1997ef about 11 days before 
maximum brightness(^icfc line) is compared with two synthetic spectra calculated with 
the ML MONTE CARLO code. The thin solid line and the dotted line are for models having 
kinetic energies of 8 x 10 51 and 1.75 x 10 52 erg, respectively. (From Q) 



well for the early circumstellar-dominated spectra (Fig. 13) . The inferences from 
spectroscopy are that the progenitor of SN 1998S underwent several distinct 
mass loss episodes during the decades and centuries before explosion, and that 
the SN ejecta or the circumstellar matter (or both) were significantly asym- 
metric. A detailed study of the nebular spectra of the Type Iln SN 1995N has 
revealed the existence of three distinct kinematic components (Fig. 14 shows 
nebular phase synthetic spectra calculated for the ~ 5000 km s _1 intermediate- 
velocity component that corresponds to unshocked SN ejecta) and produced the 
suggestion that SNe Iln are produced by red supergiants that collapse during 
their superwind phases Q (see also Jl9|). 

Synthetic spectra have not yet been compared to spectra of hyperenergetic 
SNe Iln such as SNe 1997cy g|,[l24) and SN 1999E §7|, both of which may have 
been associated with gamma ray bursts. These events may prove to be related 
to SN 1988Z, which has been modeled in terms of CSI involving either dense 
clumps or a dense equatorial wind p3 . 
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Fig. 12. The spectrum of the Type II SN 1999em four days before maximum brightness 
(thin line) is compared to a PHOENIX synthetic spectrum (thick line) with Z = Zq/IOO. 
All lines in the synthetic spectrum are produced by hydrogen and helium. The inset 
shows a shoulder in both the observed and calculated H/3 profiles. (From [^]) 

7 Prospects 

Observationally, in addition to a large increase in the sheer number of optical 
spectra, we can look forward to more and better infrared and nebular spectra, 
spectra obtained shortly after explosion, and spectropolarimetry. 

Fast spherically symmetric (ID) photospheric-phase codes such as SYNDW and 
ML MONTE CARLO will continue to be valuable for making rapid interpretations 
of new flux spectra and for certain comparative studies. Monte Carlo ID NLTE 
codes may be constructed |8q| . As computational resources permit, detailed ID 
codes such as PHOENIX will be used to generate large grids of models and spectra 
to expedite the comparison with observed spectra. PHOENIX will be modified to 
take into account the radiative coupling between circumstellar matter and SN 
ejecta, and as forbidden lines are added to the line list PHOENIX also will become 
a powerful ID nebular-phase code. 

The major challenge for the future is to develop codes for calculating spec- 
tra for arbitrary geometries (3D) of both the SN ejecta and the circumstellar 
matter (and learning how to use such codes effectively) . For flux spectra, the al- 
ready existing CLUMPYSYN code will be followed in the near future by 3D Monte 
Carlo codes and, eventually, by 3D PHOENIX-level codes. Spectropolarization 
calculations for 3D configurations are certain to emerge as a major focus of SN 
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Fig. 13. The UV (black line) and optical {blue line) spectra of the Type Iln SN 1998S 
about four days before maximum brightness are compared with a PHOENIX synthetic 
spectrum (red line) calculated for a circumstellar shell having an r~ 2 density distribu- 
tion and a constant expansion velocity of 1000 km s -1 . (From LenOla) 



spectroscopy; Monte Carlo codes for calculating 3D polarization spectra also will 
become available in the near future. 

We are grateful to Doug Leonard for providing an unpublished figure, and to 
him and Lifan Wang for helpful comments. This work was supported by NASA, 
NSF, and the Physics Department of New Mexico Tech. 
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Rest wavelength (A) 

Fig. 14. The spectrum of the Type Iln SN 1995N nearly five years after maximum 
brightness is compared to synthetic spectra calculated with (upper) and without (lower) 
the influence of Lya fluorescence on the Fe II lines taken into account. (From J42|] ) 
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